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Abstract. Using our interferometric angular diameter measurements of seven classical Cepheids reported in 
Kervella et al. (2003 Paper I), complemented by previously existing measurements, we derive new calibrations 
of the Cepheids Period-Radius (P-R) and Period-Luminosity (P-L) relations. We obtain a P-R relation of 
log R = [0.767 ± 0.009] log P + [1.091 ± 0.011], only 1 a away from the relation obtained by Gieren et al. 
We therefore confirm their P-R relation at a level of A(log7?) = ±0.02. We also derive an original calibration of 
the P-L relation, assuming the slopes derived by Gieren et al. (19981 from LMC Cepheids, ax = —3.267 ± 0.042 
and av = -2.769 ± 0.073. With a P-L relation of the form Mx = ctx (logP - 1) + /3a, we obtain logP = 1 
reference points of /3a- ~ —5.904 ± 0.063 and /?v = —4.209 ± 0.075. Our calibration in the V band is statistically 
identical to the geometrical result of Lanoix et al. 11999t . 
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1. Introduction 

The Period-Luminosity (P-L) relation of the Cepheids is 
the basis of the extragalactic distance scale, but its cal- 
ibration is still uncertain at a AAI — ±0.10 mag level. 
Moreover, it is not excluded that a significant bias of the 
same order of magnitude affects our current calibration of 
this relation. On the other hand, the Period-Radius re- 
lation (P-R) is an important constraint to the Cepheid 
models (see e.g. Alibert et al. 11999(1 . 

Traditionally, there has been two ways to calibrate the 
P-L relation. For Cepheids in clusters one can use main 
sequence fitting, assuming that the main sequence is sim- 
ilar to that of the Pleiades. This method has been ques- 
tioned however, following the release of Hipparcos data 
(e.g. Pinsonneault et al. 119981 but see also Pan et al. 120041 
and Robichon et al. I1999|l . Another route to the P-L re- 
lation is the Baade-Wesselink (BW) method where one 
combines photometry and radial velocity data to obtain 
the distance and radius of a Cepheid. Recent applications 
of the BW method to individual stars can be found for in- 
stance in Taylor et al. lf1W7|l and Taylor & Booth (1998;), 
while the calibration of the P-R and P-L relations using 
BW distances and radii is demonstrated in Gieren, Fouque 
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& Gomez IjlMSl hereafter GFG98). A requirement of this 
method is a very accurate measurement of the Cepheid's 
effective temperature at all observed phases, in order to 
determine the angular diameter. Interferometry allows us 
to bypass this step and its associated uncertainties by 
measuring directly the variation of angular diameter dur- 
ing the pulsation cycle. As shown by Kervella et al. H2003L 
hereafter Paper I) and Lane et al. ( 2002 ) , the latest gener- 
ation of long baseline visible and infrared interferometers 
have the potential to provide precise distances to Cepheids 
up to about 1 kpc, using the interferometric BW method 
(see Sect. 121. 

The main goal of the present paper is to explore the 
application of this technique to the calibration of the P- 
R and P-L relations, and to verify that no large bias is 
present in the previously published calibrations of these 
important relations. Our sample is currently too limited 
to allow a robust determination of the P-L relation, de- 
fined as A/a = OA (log P — 1) ± /3a, that would include 
both the slope a\ and the logP = 1 reference point /3a- 
However, if we suppose that the slope is known a priori 
from the literature, we can derive a precise calibration of 
/3a- In Sect. 121 we present our determination of the P-R re- 
lation using new angular diameter values from Paper I, as 
well as previously published interferometric and trigono- 
metric parallax measurements. Sect. 0] is dedicated to the 
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calibration of the P-L relation reference points f3\ in the 
K and V bands. The consequences on the LMC distance 
are briefly discussed in Sect. 14.51 

2. Cepheid distances by interferometry 

We have obtained angular diameter measurements for 
seven Cepheids with the VLT interferometer (Kervella 
et al. 2003, Paper 1). These if -band measurements were 
made with the VINCI instrument (Kervella et al. I2()()3|l 
fed by two 0.35m siderostats. Several basehnes were used, 
ranging from 60 m to 140 m. Our measurements, described 
in detail in Paper I, have a typical precision of 1 to 3 %. 
This is good enough to actually resolve the pulsation of 
several Cepheids; in other words we can follow the change 
in angular diameter. We have combined these measure- 
ments with radial velocity data and derived a radius and 
distance for four Cepheids of our sample. For the remain- 
ing three stars, we were able to derive their mean angular 
diameters, but the pulsation remained below our detec- 
tion threshold. This sample was completed by previously 
published measurements obtained with other instruments. 

In the present work, we have retained the limb dark- 
ened (LD) angular diameters 0ld provided by each author. 
Marengo et al. (|2002,, 2,003| have shown that the LD prop- 
erties of Cepheids can be different from those of stable 
stars, in particular at visible wavelengths. For the mea- 
surements obtained using the GET (Mourard et al. ll997|l 
and NPOI (Nordgren et al. 2000), the LD correction is 
relatively large {k — 6'ld/6'ud — 1.05), and this could be 
the source of a bias at a level of a 1 to 2% (Marengo et 
al. l2004|l . However, in the infrared, the correction is much 
smaller {k ~ 1-02), and the error on its absolute value is 
expected to be significantly below 1%. The majority of the 
Cepheid interferometric measurements was obtained in 
the H and K bands (FLUOR/IOTA, PTI, VLTI/VINCI), 
and we believe that the potential bias introduced on our 
fits is significantly smaller than their stated error bars. 
The final answer about the question of the limb darkening 
of Cepheids will come from direct interferometric obser- 
vations. The direct measurement of the limb darkening of 
nearby Cepheids will soon be possible with the AMBER 
instrument (Petrov et al. l^nnHjl of the VLTI. 

The radial velocity data were taken from Bersier 
(2002). They have been obtained with the CORAVEL 
spectrograph (Baranne, Mayor & Poncet, 1979). This in- 
strument performs a cross-correlation of the blue part of a 
star's spectrum (3600 — 5200 A) with the spectrum of a red 
giant. A gaussian function is then fitted to the resulting 
cross-correlation function, yielding the radial velocity. 

In Paper I, we have applied three distinct methods (or- 
ders 0, 1 and 2) to derive the distances d to seven Galactic 
Cepheids from interferometric angular diameter measure- 
ments. Not all three methods can be used to derive the 
distance for every star, depending on the level of com- 
pleteness and precision of the available angular diameter 
measurements : 



— Order 0: constant diameter model. 

This is the most basic method, used when the pulsation 
of the star is not detected. The average linear diameter 
D of the star is supposed to be constant and known a 
priori, e.g. from a previously published P-R relations 
(such as the relation derived by GFG98). Knowing the 
linear and angular radii, the only remaining variable 
to fit is the distance d. 

— Order 1: variable diameter model. 

We still consider that the average linear diameter of the 
star is known a priori, but we include in our angular 
diameter model the radius variation curve derived from 
the integration of the radial velocity of the star. This 
method is well suited when the intrinsic accuracy of the 
angular diameter measurements is too low to measure 
precisely the pulsation amplitude. The distance d is 
the only free parameter for the fit. 

— Order 2: interferometric BW method. 

The interferometric variant of the BW method 
(Davis ,19791 Sasselov et al. 11994(1 combines the an- 
gular amplitude of the pulsation measured by inter- 
ferometry and the linear displacement of the stellar 
photosphere deduced from the integration of the ra- 
dial velocity curve to retrieve the distance of the star 
geometrically. This method is also called "parallax of 
the pulsation". In the fitting process, the radius curve 
is matched to the observed angular diameter curve, us- 
ing both the distance and linear diameter as variables. 
Apart from direct trigonometric parallax, this method 
is the most direct way of measuring the distance of a 
Cepheid. It requires a high precision angular diameter 
curve and a good phase coverage. 

The order 0/1 methods, on one hand, and 2 on the 
other hand, are fundamentally different in their assump- 
tions, and the distance estimates are affected by different 
kinds of errors. While the order 2 method errors are due 
to the interferometric measurement uncertainties (mostly 
statistical), the order 0/1 distances carry the systematic 
error bars of the assumed P-R relation. As they are fully 
correlated for all stars in the sample, they cannot be aver- 
aged over the sample. In particular, the order 0/1 diame- 
ters cannot be used to calibrate the P-R relation, as they 
assume this relation to be known a priori. 

Due to its stringent requirements in terms of precision, 
the interferometric BW method (order 2) was applied suc- 
cessfully up to now on five Cepheids only: (Car (Paper I), 
/5Dor (Paper I), jyAql (Paper I; Lane et al. .2002J . WSgr 
(Paper I) and ( Gem (Lane et al. l2002|l . However, it is ex- 
pected that many more stars will be measurable with the 
required precision in the near future (see Sect.EJ. 

3. Period-Radius relation 

3.1. Method 

The Period-Radius relation (P-R) of the Cepheids takes 
the form of the linear expression: 

logR^a\ogP + b (1) 
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Fig. 1. Period-Radius diagram deduced from the inter- 
ferometric observations of Cepheids listed in Table ^ 
The thin dashed hne represents the best-fit P-R 
relation assuming the slope of GFG98, logi? = 
0.750 [±0.024] logP -F 1.105 [±0.017 ± 0.023]. The solid 
line is the best-fit relation allowing both the slope and 
zero point to vary, logi? = 0.767 [±0.009] logP + 
1.091 [±0.011]. 



In order to calibrate this relation, we need to estimate 
directly the linear radii of a set of Cepheids. We have ap- 
plied two methods to determine the radii of the Cepheids 
of our sample: the interferometric BW method, and a com- 
bination of the average angular diameter and trigonomet- 
ric parallax. While the first provides directly the average 
linear radius and distance, we need to use trigonometric 
parallaxes to derive the radii of the Cepheids for which 
the pulsation is not detected. We applied the Hipparcos 
parallaxes (Perryman et al. ll997|l to all the order 0/1 mea- 
surements, except (S Cep, for which we considered the re- 
cent measurement by Benedict et al. (|2002|l . Table Q] lists 
the Cepheid linear radii that we obtain. 

We can use the results from both order 0/1 and 2 meth- 
ods at the same time, as the obtained linear radii obtained 
in this way are fully independent from each other. On 
one hand (BW method), we obtain them considering the 
amplitude of the pulsation and the radial velocity curve, 
while on the other hand, they are derived from the aver- 
age angular diameter and the trigonometric parallax. As 
the amplitude of the pulsation and the average diameter 
values are distinct observables, these two methods can be 
used simultaneously in the fit. 



Table 2. Period-Radius relations, assuming an expres- 
sion of the form logi? — a \ogP + b. For the fitting of 
b alone, the slope has been assumed as known a priori 
from GFG98. In this case, its error bar translates to a 
systematic uncertainty on the b value derived from the fit 
(given in brackets). References: (1) GFG98, (2) Turner & 
Burke PUH^ . (3) This work. 



Ref. 


Fit 


a ± (Tstat 


b ± (Tstat [±0"syst] 


(1) 




0.750 ± 0.024 


1.075 ± 0.007 


(2) 




0.747 ± 0.028 


1.071 ±0.025 


(3) 


b only 




1.105 ±0.017 [±0.023] 


(3) 


a, b 


0.767 ±0.009 


1.091 ±0.011 



tematic errors). As a comparison, GFG98 have obtained 
a value of b = 1.075 ± 0.007, only -1.6 a away from our 
result. The relations found by Turner & Burke H2002|l and 
Laney & Stobie p995 P are very similar to GFG98, and 
are also compatible with our calibration within their error 
bars. 

Fitting simultaneously both the slope and the zero 
point to our data set, we obtain a = 0.767 ± 0.009 and 
b = 1.091 ± 0.011. These values are only Aa = +0.7 a 
and Ah = ±1.2 tr away from the GFG98 calibration. 
Considering the limited size of our sample, the agreement 
is very satisfactory. On the other hand, the slopes derived 
by Ripepi et al. (119971 and Krockenberger, Sasselov & 
Noyes (1997), both around 0.60, seem to be significantly 
too shallow. 



4. Period-Luminosity relation 

4.1. Distance estimates 

For the order and 1 methods (Paper 1) , we used an a pri- 
ori P-R relation (from GFG98) to predict the true linear 
diameter of the Cepheids of our sample. This relation re- 
lies on the measurement of the photometric flux, effective 
temperature (classical BW method) and radial velocity. 
The apparent magnitude also intervenes in the computa- 
tion of the absolute magnitude, and therefore we cannot 
use these distance estimates to calibrate the P-L relation 
without creating a circular reference. For this reason, we 
have considered only the distances obtained using the in- 
terferometric BW method (order 2) for our P-L relation 
calibration, complemented by the Benedict et al. H2002|l 
trigonometric parallax of S Cep. 



3.2. Calibration results 

Fig. n shows the distribution of the measured diame- 
ters on the P-R diagram, based on the values listed in 
Tabled When we choose to consider a constant slope of 
a = 0.750 ± 0.024, as found by GFG98, we derive a zero 
point of 6 = 1.105 ± 0.017 ± 0.023 (statistical and sys- 



4.2. Absolute magnitudes 

The average apparent magnitudes in V and K of S Cep 
were computed via a Fourier series fit of the data from 
Moffett & Barnes (TlgSgl and Barnes et al. for the 

K band and Barnes et al. (jl997|l for the V band. The 
sources for the other apparent magnitudes are given in 
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Table 1. Weighted averages of the interferometric mean angular diameters (?ld and of the geometric distances d to 
nearby Cepheids (bold characters). These values were used to compute the linear radii given in the last two columns. 
The individual measurements used in the averaging process are also given separately for each star. References: (1) 
Mourard et al. (f1WZ|, (2) Nordgren et al. (SOOT), (3) Lane et al. ipHII^jl . (4) Mozurkewich et al. l(TMT|l . (5) Paper I, 
(6) Benedict et al. (2^), (7) Ferryman et al. (1997). 



Star 


P(d) 


logP 


Ref. 6'ld 


6'ld (mas) 


Ref. d 


d(pc) 


R{Rq) 


logi? 




S Cap 


5.3663 


0.7297 


(1) 


1.521 ±0.010 
1.60 ±0.12 




274+i^ 


44.8li:^ 


1.651: 


-0.018 








(2) 


1.52 ±0.01 




















(6) 
(7) 










XSgr 


7.0131 


0.8459 




1.471 ±0.033 




330±^f 


52.2+^^ 


1.717; 


fU.lBl 
-0.118 








(5) 


1.471 ±0.033 




















(7) 


330l^f 








7] Aql 


7.1768 


0.8559 


(2) 
(3) 
(5) 


1.791 ±0.022 

1.69 ±0.04 
1.793 ± 0.070 
1.839 ±0.028 


(3) 
(5) 
(7) 


308+24 

320tg 
276tl 


59-3+4.a 


1.773; 


fU.U37 
-0.035 


W Sgr 


7.5949 


0.8805 


(5) 


1.312 ±0.029 

1.312 ±0.029 


(5) 
(7) 


on-q+216 
"J ' ^—130 

637l^f, 


56.4l?^ 


1.751; 


fU.184 
-0.146 


P Dor 


9.8424 


0.9931 


(5) 


1.884 ±0.024 

1.884 ±0.024 


(5) 
(7) 


323+42 

3451^^^ 

31811;; 


65.4+- 


1.816; 


-0.061 


C Gem 


10.1501 


1.0065 


(2) 
(3) 
(4) 

(5) 


1.688 ±0.022 
1.55 ±0.09 
1.675 ± 0.029 
1.73 ±0.05 
1.747 ±0.061 


(3) 
(7) 


362+j^^ 

362-tll 

358lif 


65.6+-^ 


1.817; 


-0.044 


Y Oph 


17.1269 


1.2337 




1.438 ±0.051 




'-360 


136+53" 


2.132; 


fU.63i 
-0.231 








(5) 


1.438 ±0.051 




















(7) 


077+2100 
8' '-360 








£ Car 


35.5513 


1.5509 




2.988 ±0.012 




597+^^ 


191.2+;:^ 


2.281; 


fU.UlT 

-0.014 








(5) 


2.988 ±0.012 


(5) 
(7) 


603+?^ 
463i-^ 





Paper I (Table 1). Following Fouque et al. pM| . the 
extinction Ax has been computed using the relations: 

Ax = Rx Eb-v (2) 

i?v = 3.07 ±0.28(5 -y)o± 0.04 £;B_y (3) 

i?^ = i?y /II ~ 0.279 (4) 

The resulting extinction values are listed in Table O and 
the final absolute magnitudes Mx of the Cepheids of our 
sample are listed in Table El 

4.3. Calibration of tine P-L relation 

We have considered for our fit the P-L slope measured 
on LMC Cepheids. This is a reasonable assumption, as 
it can be measured precisely on the Magellanic Clouds 
Cepheids, and in addition our sample is currently too lim- 
ited to derive both the slope and the logP = 1 reference 
point simultaneously. 



Table 3. Apparent magnitudes and extinctions in the 
K and V bands for the Cepheid whose distances have 
been measured directly by interferometry. {B — V)q is the 
mean {B — V) index as reported in the online database 
by Fernie et al. H1995(l . The Eb-v values were taken from 
Fernie H1990|l . The extinctions in the K and V bands are 
given respectively in the " Ak" and "Ay" columns, in mag- 



nitudes. 


Star 


(B - V)o 


Eb-v 


ruK 


Ak 


mv 


Av 


SCep 


0.66 


0.09 


2.31 


0.03 


3.99 


0.30 


?7 Aql 


0.79 


0.15 


1.97 


0.04 


3.94 


0.49 


WSgr 


0.75 


0.11 


2.82 


0.03 


4.70 


0.36 


/3Dor 


0.81 


0.04 


1.96 


0.01 


3.73 


0.15 


CGem 


0.80 


0.02 


2.11 


0.01 


3.93 


0.06 


^Car 


1.30 


0.17 


1.09 


0.05 


3.77 


0.58 



Recently, Fouque et al. H2003|l have revised the P-L 
slopes derived from the large 0GLE2 survey (Udalski et 
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Table 4. Absolute magnitudes of Cepheids measured 
exclusively using the interferometric Baade-Wesselink 
method, except for 6 Cep, whose parallax was taken from 
Benedict et al. The same error bars apply to the 

K and V band absolute magnitudes. The Cepheid periods 
are listed in Tabled References: (1) Lane et al. I|2U(J2|I . 
(2) Benedict et al. (PUH^ . (3) Paper I. 



Star Ref. d ±a Mk My ±(t 



(5 Cep 


(2) 


273 


+ 12 
-11 


-4.90 


-3.49 


+0.09 
-0.09 


Aql 


(1) 


320 


+32 
-32 


-5.60 


-4.08 


+0.23 
-0.21 


J? Aql 


(3) 


276 


+55 
-38 


-5.28 


-3.76 


+0.32 
-0.39 


WSgr 


(3) 


379 


+216 
-130 


-5.10 


-3.56 


+0.91 
-0.98 


/3Dor 


(3) 


345 


+ 175 
-80 


-5.74 


-4.10 


+0.57 
-0.89 


CGem 


(1) 


362 


+38 
-38 


-5.69 


-3.92 


+0.24 
-0.22 


£Car 


(3) 


603 


+24 
-19 


-7.86 


-5.72 


+0.07 
-0.08 



al. lTM^ . and obtain values of ay = -2.774 ± 0.042 and 
UK = —3.215 ± 0.037. These values are consistent within 
their error bars with LPG99 {ay = -2.77 ±0.08), GFG98 
(ay = -2.769±0.073, uk = -3.267±0.042) and Sasselov 
et al. 1)19971 ay = — 2.78±0.16). Considering this consen- 
sus, we have chosen to use the slope from GFG98 to keep 
the consistence with the P-R relation assumed in Paper I. 

Tables and El report the results of our calibrations 
of the P-L relations, and the positions of the Cepheids on 
the P-L diagram are shown on Fig.|21and Fig.O The final 
log P = 1 reference points are given in bold characters in 
Tables |5l and Our calibrations differ from GFG98 by 
ABk — -1-0.20 mag in the K band, and Aby = -1-0.14 mag 
in V, corresponding to -1-3.0 and -1-1.8 ct, respectively. The 
sample is dominated by the high precision £ Car and 6 Cep 
measurements. When these two stars are removed from 
the fit, the difference with GFG98 is shghtly increased, up 
to -1-0.25 and -1-0.30 mag, though the distance in a units is 
reduced (-1-1.3 and -1-1.5). From this agreement, .^Car and 
S Cep do not appear to be systematically different from 
the other Cepheids of our sample. 

It is difficult to conclude firmly to a significant dis- 
crepancy between GFG98 and our results, as our sample 
is currently too limited to exclude a small-statistics bias. 
However, if we assume an intrinsic dispersion of the P-L 
relation ctpl — 0.1 mag, as suggested by GFG98, then our 
results point toward a slight underestimation of the ab- 
solute magnitudes of Cepheids by these authors. On the 
other hand, we obtain precisely the same log P — 1 refer- 
ence point value in V as Lanoix et al. (19211 using paral- 
laxes from HiPPARCOS). The excellent agreement between 
these two fully independent, geometrical calibrations of 
the P-L relation is remarkable. 

4.4. P-L relation slopes in the Galaxy and in the LMC 

The question of the difference in slope between the 
Galactic and LMC Cepheid P-L relations has recently 
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Fig. 2. Period-Luminosity diagram in the K band using 
only interferometric BW distances and the <5 Cep parallax 
listed in Table0| The solid line represents the best-fit P-L 
relation using the slope derived by GFG98 (classical least- 
squares fit: the individual measurements are weighted by 
the inverse of their variance) . 

Table 5. Period-Luminosity relation intercept (ix for a 
10 days period Cepheid (log P = 1), in the K band. We as- 
sume an expression of the form Mk = olk (log P—X)+ (3 k ■ 
The slope value is taken from GFG98 {aK = -3.267 ± 
0.042). The systematic error is given in brackets, and cor- 
responds to the uncertainty on the GFG98 slope. 



Ref. 




icTstat 


±0"syst 


GFG98 


-5.701 


±0.025 




This work, all stars 


-5.904 


±0.063 


±0.005 


Without 5 Cep and £ Car 


-5.956 


±0.191 


±0.006 



been discussed by Fouque et al. H2()()3|l and Tammann et 
al. 1)2003(1 . These authors conclude that the Galactic slopes 
are significantly steeper than their LMC counterparts. 
For example, Tamman et al. f20035 obtain ay [Gal] — 
-3.14 ±0.10, while Fouque et al. (j2003iil derive ay [Gal] = 
-3.06 ± 0.11 and ay [LMC] = -2.774 ± 0.042. 

However, our fit is largely insensitive to the precise 
value assumed for the P-L relation slope. Considering 
the steeper Tammann et al. Ij2003|) slope, we obtain a 
best fit logP = 1 absolute magnitude of f3y = —4.211 ± 
0.075 ± 0.001, identical to the calibration obtained using 
the GFG98 slope. The small systematic error bar that we 
obtain on /3y (corresponding to the ±0.10 error on ay) 
shows the weakness of the correlation between a and P in 
our fit. However, the reduced of the fit is significantly 
larger with this steeper slope (Xred ~ 1-25) than with the 
LMC slope from GFG98 (xL = 0-53). 
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Fig. 3. Period-Luminosity diagram in the V band (slope 
from GFG98). 



Table 6. Period-Luminosity relation intercept /3y 
(logP = 1) in the V band, derived using the GFG98 slope 
(ay = -2.769 ± 0.073). 





/3v 


io"stat 


±C'"syst 


GFG98 


-4.063 


±0.034 




LPG99 


-4.21 


±0.05 




This work, all stars 


4.209 


±0.075 


±0.001 


Without 5 Cep and I Car 


-4.358 


±0.197 


±0.010 



4.5. The distance to the LMC 

The apparent magnitudes in V and if of a 10 days pe- 
riod Cepheid in the Large Magellanic Cloud (LMC) de- 
rived by Fouque et al. (2003) from the OGLE Cepheids 
are ZP^ = 12.806 ± 0.026 and ZPy = 14.453 ± 0.029. 
These authors assumed in their computation a constant 
reddening of E{B ~ V) ^ 0.10 for all the LMC Cepheids 
they have used (more than 600). Our calibrations of the 
Galactic Cepheids P-L relations in K and V thus im- 
plies LMC distance moduli of /.ik ~ 18.71 ± 0.07 and 
yUy = 18.66 ± 0.08, respectively. 

From a large number of photometric measurements of 
LMC and SMC Cepheids obtained in the framework of 
the EROS programme, Sasselov et al. H1997|l have shown 
that a Sfi correction has to be applied to the LMC dis- 
tance modulus to account for the difference in metallicity 
between the LMC and the Galactic Cepheids. They have 
determined empirically a value of: 

= Mtruo - Mobscrvod = -0.14 ± 0.06 (5) 

We would like to point out that this correction has been 
questioned by Udalski et al. H2UU1|I . based on Cepheid ob- 
servations in a low metallicity galaxy (IC1613), and its 
amplitude is still under discussion (Fouque et al. I2003|l . 

Averaging our K and V band zero point values (with- 
out reducing the uncertainty, that is systematic in na- 



ture), we obtain a final LMC distance modulus of /iq = 
18.55 ± 0.10. This value is only +0.8 a away from the 
Mo = 18.46 ± 0.06 value obtained by GFG98, and -Icr 
from the fj,o = 18.70 ± 0.10 value derived of Feast & 
Catchpole l|1997|l . It is statistically identical to the LMC 
distance used by Freedman et al. H2001|l for the HST Key 
Project, fiQ = 18.50 ± 0.10. Alternatively, if we consider 
the smaller metallicity correction of 6^ — 0.06 ± 0.06 
proposed by GFG98, we obtain a distance modulus of 
fio = 18.63 ±0.10. 

5. Conclusion and perspectives 

We have confirmed in this paper the P-R relation of 
GFG98 and Turner & Burke piH^ . to a precision of 
A(logi?) = ±0.02. We also derived an original calibra- 
tion of the P-L relations in K and V, assuming the slopes 
from GFG98 that were established using LMC Cepheids. 
Our P-L relation calibration yields a distance modulus of 
fj,o = 18.55 ± 0.10 for the LMC, that is statistically iden- 
tical to the value used by Freedman et al. (2001) for the 
HST Key Project. We would like to emphasize that this re- 
sult, though encouraging, is based on six stars only (seven 
measurements, dominated by two stars), and our sample 
needs to be extended in order to exclude a small-number 
statistics bias. In this sense, the P-L calibration presented 
here should be considered as an intermediate step toward 
a final and robust determination of this important relation 
by interferometry. 

While our results are very encouraging, the calibra- 
tion of the PR and PL relations as described here may 
still be affected by small systematic errors. In particular 
the method relies on the fact that the displacements mea- 
sured through interferometry and through spectroscopy 
(integration of the radial velocity curve) are in different 
units (milli-arcseconds and kilometers respectively) but 
are the same physical quantity. This may not be the case. 
The regions of a Cepheid's atmosphere where the lines 
are formed do not necessarily move homologously with 
the region where the if -band continuum is formed. This 
means that the two diameter curves may not have exactly 
the same amplitude; there could even be a phase shift be- 
tween them. As discussed in Sect. |21 the limb darkening 
could also play a role at a level of ~ 1%. A full exploration 
of these effects is far beyond the scope of this paper. We 
can nevertheless put an upper bound on the systematic er- 
ror that could result from this mismatch. Our PL relation 
can be compared to that derived from Cepheids in open 
clusters, whose distances are obtained via main sequence 
fitting. The two distance scales are in excellent agreement 
(Gieren & Fouque 1993; Turner & Burke, 2002). These 
distances are consistent with a Pleiades distance modulus 
of 5.56; if anything they are slightly larger. 

The availability of 1.8 m Auxiliary Telescopes (Koehler 
et al. nU(W i on the VLTI platform in 2004, to replace 
the current 0.35 m Test Siderostats, will allow to observe 
many Cepheids with a precision at least as good as the 
observations of £ Car reported in Paper I (angular diame- 
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ters accurate to 1%). In addition, the AMBER instrument 
(Petrov et al. I2000|l will extend the VLTI capabilities to- 
ward shorter wavelengths (J and H bands), thus provid- 
ing higher spatial resolution than VINCI {K band). The 
combination of these two improvements will extend signif- 
icantly the accessible sample of Cepheids, and we expect 
that the distances to more than 30 Cepheids will be mea- 
surable with a precision better than ±5 %. This will pro- 
vide a high precision calibration of both the logP = 1 
reference point (down to ±0.01 mag) and the slope of 
the Galactic Cepheid P-L. As the galaxies hosting the 
Cepheids used in the Key Project are close to solar metal- 
licity in average (Feast ''iOOlJI, this Galactic calibration will 
allow to bypass the LMC step in the extragalactic dis- 
tance scale. Its attached uncertainty of ±0.06 due to the 
metallicity correction of the LMC Cepheids will therefore 
become irrelevant for the measurement of Hq. 
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